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We study the prototype 5 d pyrochlore iridate Y 2 lr 207 from first principles using the local density 
approximation and dynamical mean-field theory (LDA+DMFT). We map out the phase diagram in 
the space of temperature, onsite Coulomb repulsion U, and filling. Consistent with experiments, we 
find that an all-in-all-out ordered insulating phase is stable for realistic values of U. The trigonal 
crystal field enhances the hybridization between the j e ff =1/2 and j e ff = 3/2 states, and strong 
inter-band correlations are induced by the Coulomb interaction, which indicates that a three-band 
description is important. We demonstrate a substantial band narrowing in the paramagnetic metallic 
phase and non-Fermi liquid behavior in the electron/hole doped system originating from long-lived 
quasi-spin moments induced by nearly flat bands. 

PACS numbers: 71.15.Mb,71.27.+a,71.30.+h 


The competition and cooperation between spin-orbit 
coupling (SOC) and electron correlations induces novel 
phenomena in 4 d and 5 d transition metal oxides such as 
spin-orbit-assisted Mott insulators, topological phases, 
and spin liquids [I]. The pyrochlore iridates 
(A= Pr, Nd, Y, etc.) are an ideal system to study these 
phenomena because their magnetic and electronic states 
can be tuned by chemical substitution, pressure, and 
temperature (T). Furthermore, intriguing phenomena 
such as correlated topological phases have been theoret¬ 
ically predicted on their geometrically frustrated crystal 
structure pQ. 

In 2001, it was reported that these compounds show 
a crossover from metal to insulator with decreasing A 3+ 
ionic radii at high T [2] and a magnetic anomaly was 
found at low T for small A 3+ ionic radii [3]. For the 
metallic compound A=Pr, experiments revealed spin- 
liquid behavior [4j [5] and an unconventional anomalous 
Hall effect [6]. On the other hand, the Ir magnetic or¬ 
dering has not been determined for a decade due to the 
strong neutron absorption by Ir and large magnetic con¬ 
tributions from rare-earth / moments on A 3+ . The mag¬ 
netic order has only recently been identified as a non- 
collinear all-in-all-out order [see Fig. 0a)] CHj. 

Among the insulating compounds, Y2R2O7 has the 
highest magnetic transition temperature and no / mo¬ 
ments. This makes this compound a prototype sys¬ 
tem for studying strong electron correlations among 5 d 
electrons. A pioneering local density approximation 
(LDA)+E/ study for this compound showed that the all- 
in-all-out order is indeed stable at large on-site repulsion 
u HO]. It also proposed a topological Weyl semimetal 
as the ground state of some compounds in this series. 
This stimulated further theoretical studies on the topo¬ 
logical nature and unconventional quantum criticality of 
5 d electrons on the pyrochlore lattice EDUMg. 


In pyrochlore iridates, the Ir atoms form a frustrated 
pyrochlore lattice, a corner-sharing network of tetrahedra 
[see Fig. [lja) and Ref. 20 . The so-called j e s=l/2 pic¬ 
ture was originally proposed for the insulating quasi-2D 
compound Sr2lr04 with the same electron configuration 
5 d 5 (21j [22]. The SOC splits the 1 2 g manifold into a 
fully occupied j e ff=3/2 quartet and a half-filled j e fr=l/2 
doublet = S — L is the effective total angular mo¬ 
mentum, S and L the spin and orbital momenta). This 
picture was subsequently confirmed by LDA+dynamical 
mean-field theory (DMFT) studies [231425] . The j e ff=l/2 
and 3/2 manifolds hybridize either under a tetragonal 
crystal field (CF), as in the case of Sr 2 Ir04, or under a 
trigonal CF because they do not commute with j 2 ff . For 
pyrochlore iridates, a recent quantum chemistry calcu¬ 
lation found that the trigonal CF is comparable to the 
SOC [26]. Moreover, the t 2 g band width in the LDA band 
structure (~ 2 eV W) is even smaller than the typical 
value of U for Ir 4+ (~ 3 eV [28]). These effects will en¬ 
hance the hybridization between the two j e ^ manifolds 
beyond the single-particle level, raising questions about 
the validity of the j e ft= 1/2 single-band picture. Another 
interesting issue is the role of the geometrical frustra¬ 
tion, which suppresses long-range magnetic ordering and 
induces a novel type of magnetism [20]. To address these 
questions, it is important to perform a first-principles 
study of the three-band system which takes into account 
SOC, electron correlations, itinerancy, and the crystal 
structure. 

In this Letter, we present state-of-the-art relativistic 
LDA+DMFT calculations for the prototypical compound 
Y2lr2 07. We map out the phase diagram in T and on¬ 
site Coulomb repulsion V and investigate the properties 
of the correlated 5 d electrons in this family. First, we con¬ 
struct maximally localized Wannier functions for the t 2 g 
manifold [29] using an LDA exchange-correlation func- 
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FIG. 1. (color online), (a) Pyrochlore lattice formed by Ir 
atoms with arrows representing spin moments in the all-in- 
all-out magnetic structure, (b) fee unit cell with the local co¬ 
ordinate axes and the energy diagram under SOC (£) and the 
trigonal crystal field (A tr i). (c) LDA band structure together 
with the density of states projected on the j e ff basis. The bro¬ 
ken line shows the total density of states. The “j e ff=l/2” and 
“jeff=3/2” bands consists mainly of the j e s= 1/2 and j e s= 3/2 
manifolds, but are substantially hybridized under the trigonal 
CF (see the text). 


tional |30| I3T ] . We use the code QMAS (Quantum MA- 
terials Simulator) [3T, which is based on the projector 
augmented wave method [33], and the two-component 
formalism [34j [35]. The experimental crystal structure 
at 290 K is taken from Ref. [36, In our DMFT calcula¬ 
tions, electron correlation effects are taken into account 
by introducing the Slater-Kanamori interaction 
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in the standard parameterization U aaaa = U , U a p a p = 
U - 2J h , U a pp a = U aa pp = Jr (oi % /?), with a (/?) 
and a (a') being orbital and spin indices, respectively. 
U and Jr are the on-site repulsion and the Hund’s cou¬ 
pling, respectively. We choose J-r/U = 0.1, which is mo¬ 
tivated by a first-principles estimate for the related com¬ 
pound Na 2 lr 03 (17=2.72 eV, Jp = 0.23 eV) [28]. Within 
DMFT, one has to solve a three-orbital quantum impu¬ 
rity problem with off-diagonal and complex hybridization 
functions. We employ a numerically exact continuous¬ 
time quantum Monte Carlo impurity solver based on the 
hybridization expansion [32 EE]. in previous studies, the 
quantum impurity models for 5 d electrons have been sim¬ 
plified to avoid a severe sign problem, e.g., by omitting 
off-diagonal hybridization functions and some interaction 
terms in the j e ff basis [39 J. Since pyrochlore iridates have 
large inter-band hybridizations, we solve our impurity 
problem without such approximations. Another advan¬ 
tage of this treatment is that we do not necessarily need 
to assume the quantization axes of spin and orbital. The 
sign problem is reduced by rotating the single-particle 
basis of the hybridization function m- Please refer to 
the Supplemental Material for technical details and some 


results on the effects of the off-diagonal hybridizations. 

Figure [TJc) shows the computed LDA band structure. 
The upper half-filled manifold, which is usually identified 
as the jeff =1/2 manifold, has an overlap with the lower 
manifold in energy space, although the bands are sepa¬ 
rated at each k point. The j e fp=l/2 manifold has four 
Kramers degenerate bands since a unit cell contains four 
Ir atoms. We constructed a tight-binding model based on 
t 2 g -orbital-like maximally localized Wannier functions. 
The SOC £ and the trigonal crystal field A tr i are es¬ 
timated to be £ = 0.40 eV and A tri = 0.23 eV ]4T] . 
These values are consistent with an estimate by a quan¬ 
tum chemistry calculation [26]. As shown in Fig. [ljb), 
the t 2 g manifold splits into three doublets under £ and 
A t ri- The wavefunction of the highest doublet £q is given 
by </>i± = — 0.977| 1/2, ±1/2) — 0.212 13/2, ±1/2) in the j eff 
basis |Jeff ? Jeff 1 ) • We denote by the effective angular 
momentum along the local [111] axis [see Fig. 0a)]- The 
11/2, d=l/2) states have about 71% reduced spin and or¬ 
bital moments compared to the ideal atomic values 1/3/ip 
and 2/3/iB, because the Wannier functions have substan¬ 
tial weight on neighboring oxygen atoms. On the other 
hand, the magnetic moments are enhanced by the hy¬ 
bridization between the j e ff =1/2 and j e ff=3/2 manifolds 
by A t ri- As a result, the doublet £q has spin and orbital 
moments of 0.490/iB and 0.598/iB* To illustrate the ef¬ 
fects of itinerancy, we plot the density of states projected 
on the j e ff basis. The contributions of 11/2, =bl/2) and 
13/2, ±3/2), which are not mixed by A tr i, have compara¬ 
ble weight near the Fermi level, which indicates that the 
inter-atomic hybridization also plays a substantial role. 

Next, we discuss the U-T phase diagram obtained by 
the DMFT calculations [Fig. [2ja)]. There is a dome¬ 
shaped all-in-all-out magnetically ordered phase at large 
U. The transition temperature T c rises up to values 
about three times higher than the experimental T® xp < 
150 K. This may be due to the neglect of spatial fluc¬ 
tuations in the DMFT approximation. We show the In¬ 
dependent spin and orbital moments along the local [111] 
axis (mm and L m) computed at several values of U in 
Fig. [2jb) . Both order parameters emerge concurrently 
with the same sign. The Ir magnetic moment is esti¬ 
mated to be 0.93 fiB at low T for U = 2.5 eV, which is 
larger than an experimental estimate of the upper bound 
(0.5 iib) 03- The ratio Lm/mm ~ 1.2 is substantially 
smaller than Lm/mm = 2 of the j e fp=1/2 state. In con¬ 
trast, a LDA±DMFT study estimated Lm/mm ~ 2.2 
for Sr 2 lr 04 [25]. On the other hand, the spectral weight 
A(uo = 0) shows a drop at T c for U = 2 eV, signaling 
a transition from a paramagnetic metal into an all-in- 
all-out ordered insulator. At higher U (> 2.5 eV), the 
spectral function is substantially suppressed even above 
T c , indicating that the system is in a Mott insulating 
state. A crossover between a paramagnetic metal and a 
paramagnetic insulator is located around U = 2.3 eV at 
high T. In Fig. [2ja), we also show the first-order Mott 
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FIG. 2. (color online), (a) U-T phase diagram at half fill¬ 
ing. There is a hrst-order transition between the magnetic 
insulator and the paramagnetic metal at low T and small U. 
The blue shaded region denotes the hysteresis region associ¬ 
ated with this transition. The error bars reflect the uncer¬ 
tainty caused by the finite number of parameter values con¬ 
sidered. The metal-insulator crossover in the high-T param¬ 
agnetic phase is shown by a broken line. The hashed region 
represents the first-order Mott transition and its hysteresis 
region in paramagnetic DMFT calculations, (b) temperature 
dependence of the angular and magnetic moment along the 
local [111] axis and the spectral weight at u — 0 for half fill¬ 
ing. The moment values of the j e fr=1/2 and 0i doublets are 
shown by dotted and broken lines, respectively (see text). 


transition line obtained by paramagnetic DMFT calcula¬ 
tions. Apparently, the metal-insulator transition in the 
magnetic DMFT phase diagram is assisted by magnetic 
ordering. The insulating compound Y 2 lr 207 , which has 
the highest T c in the family, may be located at U ~ 2.5 
eV. 

At low T and small J7, we find a first-order transition 
between the paramagnetic metallic phase and the all-in- 
all-out ordered insulating phase (both states are topologi¬ 
cally trivial). For Nd 2 lr 207 , it was reported that its mag¬ 
netic and metal-insulator transition at T c =33 K is second 
order [42 . The first-order nature may be an artifact of 
the LDA+DMFT method. A previous LDA +U study 
found a Weyl semimetallic phase on the lower -U side of 
an all-in-all-out ordered insulating phase m- In our 
DMFT phase diagram, however, the semimetallic phase 
is taken over by the insulating phase and there is a direct 
transition between the paramagnetic metal and the mag¬ 
netic insulator. This appears to be a strong correlation 
effect. 

Figure |3] shows the spectral function A(k,uo) computed 
for T = 290 K. At U = 2 eV (paramagnetic metal), 
the upper manifold near the Fermi level shows a sub¬ 
stantial band narrowing from the LDA value of ^ 1 eV 
down to approximately 0.4 eV, while the lower manifold 
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FIG. 3. (color online). Momentum-resolved spectral function 
A(k,uj) at U = 2 eV (paramagnetic metal) and U = 2.5 eV 
(all-in-all-out ordered insulator) at 290 K. The LDA band 
structure is shown by red lines. On the right we plot the k- 
integrated spectral function projected on the j e ff basis (the 
broken line is the total spectral function). 


is smeared out by correlation effects. The low-energy 
states consist mainly of the j e fp=1/2 orbitals for u) > 0. 
A similar purification of the spectral function was found 
in LDA+DMFT studies of Sr 2 lr 04 . At U = 2 eV and 
U = 2.5 eV (all-in-all-out ordered insulator), no clear 
separation is seen between the j e s=l/2 and j e ff=3/2 
manifolds in the total spectral function A(w). 

We now investigate the effect of doping this insulating 
solution. The electron filling n is changed from 4.6 to 
5.8 (n = 5 corresponds to the undoped compound). The 
all-in-all-out order vanishes rapidly upon electron/hole 
doping at n = 4.8 and n = 5.2, respectively. We plot 
Im£(za; n ) in Fig. |4ja) for n = 4.6 and n = 5.6. In a 
Fermi liquid, ImE(zu; n ) vanishes linearly with uj n at low 
frequencies (i.e, oc with a = 1) when T is sufficiently 
low. However, for a rather wide range of frequencies, 
our data show that Im Y,(iuj n ) vanishes more slowly in 
the doped insulator. We estimated the exponent a us¬ 
ing the lowest two Matsubara frequencies assuming Im 
E(icj n ) = f3uo The result is shown in Fig. ^b) for dif¬ 
ferent n. One sees a substantial reduction of a from 
1 around n = 4.6 and n = 5.8. This non-Fermi liq¬ 
uid (NFL) or bad metallic behavior persists down to the 
lowest temperature considered (77 K). 

Theoretically, it has been reported that NFL behav¬ 
ior with a reduced power-law exponent in the frequency- 
dependent self-energy can be induced by the Hund cou¬ 
pling in multi-orbital systems [431445] . In the latter case, 
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FIG. 4. (color online), (a) frequency dependence of the imag¬ 
inary part of the self-energy projected on the j eS =l/2 man¬ 
ifold. The dotted and broken lines are oc uo n and oc u;°' 5 , re¬ 
spectively. n is the number of electrons in the fo g shell (n = 5 
is the undoped case), (b) n dependence of the exponent a 
[Im Y,(iuj n ) oc oj%]. (c) T dependence of the imaginary-time 
pseudo-spin correlation at r = /3/2 (see the text for the defi¬ 
nition). (d) spectral functions computed at 116 K. (e) optical 
conductivity Q(o;). The dotted and broken lines are oc u ~ 2 
and oc cj - 0 ' 5 , respectively. 


ImH(icj n ) shows a nonzero intersect as uj n 0 close to 
half-filling, and a fractional power-law scaling appears at 
the boundary between the NFL regime and the Fermi liq¬ 
uid regime m■ This sharp crossover has been coined the 
“spin-freezing transition”, because the scattering in the 
NFL state is related to the appearance of long-lived local 
moments. To see if the NFL behavior in the pyrochlore 
iridates has a similar origin, we plot in Fig. J4^c) the 
imaginary-time correlator (s quasi (r = /3/2)5 quasi (r = 0)) 
computed for V = 2 eV, where the quasi spin is defined 
for the highest <f>i doublet as s quasi = n^ 1+ — n^ 1 _. As 
can be seen, this correlation function extrapolates to a 
finite value as T 0 for n = 4.6 and n = 5.8, indicating 
the existence of long-lived quasi-spin moments. How¬ 
ever, these long-lived moments are not a consequence of 


the Hund coupling; they rather originate from properties 
of the band dispersion, i.e. crystal structure, which we 
confirmed by reproducing this NFL behavior even with 
Jh turned off (not shown). 

As illustrated in Fig. Qd), A(lj) shows two character¬ 
istic peaks at the upper and lower edges of the j e ff=1/2 
manifold at n = 5.2 (electron-doped Fermi liquid). These 
two peaks come from nearly flat quasi-particle bands 
in the paramagnetic phase [Fig. |3ja)], corresponding to 
peak structures seen in the LDA density of states [Fig. [l]. 
The NFL behavior appears when the Fermi level hits one 
of these peaks. This leads to the formation of local mo¬ 
ments and a reduction of the coherence temperature T*. 

One might wonder if the NFL state is stable against 
magnetic ordering other than the all-in-all-out order. 
We confirmed that the NFL state is stable by perform¬ 
ing DMFT calculations without assuming the all-in-all- 
out magnetic ordering down to 116 K. This robust¬ 
ness is likely due to the geometrically frustrated crys¬ 
tal structure which prevents long-range ordering. The 
NFL state may be observed by optical conductivity mea¬ 
surements (e.g., a = 0.5 implies an optical conductivity 
fi(cj) oc 1 /y/uj [43]), or by the T-dependence of the resis¬ 
tivity. Figure [4][e) shows the optical conductivity fi(u;) 
along [001] computed at U = 2.5 eV and 116 K. The in¬ 
dependence deviates substantially from the Fermi-liquid¬ 
like behavior <n -2 for the hole-doped case n = 4.6 around 
u) = 1 eV. On the other hand, the conductivity of the 
electron-doped compound n = 5.8 is closer to a <n -2 be¬ 
havior, which is likely due to the small self-energy [see 
Fig. [4][a)]. It was reported that the undoped metallic 
BfoL^CH shows a bad metallic behavior and the transport 
properties sensitively depend on the magnetic field [46]. 
That paper discussed the connection between these non¬ 
trivial properties and a sharp peak in the LDA density 
of states near the Fermi level. Although this compound 
has a wider band width and is considered to be weakly 
correlated, it would be interesting to investigate it with 
LDA+DMFT. 

In summary, we have mapped out the finite-T phase 
diagram of the 5 d pyrochlore iridate Y 2 lr 207 using a 
state-of-the-art relativistic LDA+DMFT approach. We 
showed that the spin and orbital moments are substan¬ 
tially enhanced by the trigonal crystal field and that 
the jeff=1/2 picture is not valid. We have also identi¬ 
fied a characteristic non-Fermi liquid self-energy which 
originates from long-lived quasi-spin moments induced 
by nearly flat bands. Due to experimental difficulties 
such as the lack of a natural cleavage plane, angle- 
resolved photoemission data of this family have been re¬ 
ported only for Bi 2 Ir 2 07 07]. However, the connection 
of the non-Fermi liquid behavior to nearly flat bands en¬ 
ables an analysis of the electronic structure by macro¬ 
scopic measurements. Although the hole-doped com¬ 
pounds Ca x Y 2 - cc Ir 2 07 show a bad metallic behavior for 
x > 0.3 [48]09], the nature of this metallic state remains 
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to be clarified. It would be interesting to explore the 
optical conductivity or T-dependent resistivity of these 
doped compounds to confirm the NFL behavior. A com¬ 
parative study with the osmate Cd 20 s 207 with a 5 d 3 
configuration would also be interesting [5QH53] . 
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Solving the impurity problem 

In this section we explain how we solved the impurity 
problem using the hybridization expansion continuous¬ 
time quantum Monte Carlo algorithm (CT-HYB). [37l 138] 
There are four iridium atoms in a unit cell. The effective 
action for the ra-th iridium atom is given by 

s = Sloe 

+ YY / dTdT ' A i°J<r'( T ' -T)ch(T')Cj<J' W, (2) 

ij a a' 0 

where A is the hybridization function [A^j a /(r) = 
A ja',ia( r )\- A(r) is a 6 x 6 matrix with non-zero off- 
diagonal elements, i and j are the index of the t 2 g -like 
Wannier functions, while a and a' denote spin. The Wan- 
niner functions are constructed in terms of the local cubic 
frame of an IrC>6 octahedron. We denote the inverse tem¬ 
perature by p. 

The local Hamiltonian consists of the one-body part 
Ho and the Slater-Kanamori interaction Hint- The one- 
body part is given by 

U 0 = P m (3) 

where Nk is the number of k points. The 24x24 matrix 
H(fc) is the LDA one-body Hamiltonian in the Wannier 
basis. Pm is a projector to the ra-th iridium atom. The 
Slater-Kanamori interaction Hint is defined in terms of 
the t 2 g -like Wannier functions. Its explicit form is given 
in Eq. (1) of the main text. 


To avoid a severe sign problem without further approx¬ 
imation, we transform the single-particle basis as 

da = Y J U* ba C b , (4) 

b 

dt = J2 u bacl (5) 

b 

where U a i) is a unitary matrix with a and b being the 
combined index of spin and orbital. In practice, we build 
the matrix U from the eigenvectors of 

Ko + SxjIs', (6) 

where S is a positive but very small constant and j^ 1 is 
the effective angular momentum projected on the local 
[111] axis. 

The second term in Eq. ©> reads 

r p 

S'hyb = Y / d-rdr'A ab {r' - t)<4(t')4(t), (7) 

a/3 Jo 

with 


A ab (r) = ^(C/ t ) ac A cd (r)U db . 


(8) 


cd 


Now, we expand the partition function Z as 


Z = Z\ 


bath 


oo .. 

E E 

n =0 ai,-)«n 

rfi rP 

/ dr, dr(■ • • / dr„dr^ 

Jo Jo 

e~ mioc T dan (t„)4, (t'J ■ ■ ■ d ai (ri)dL (t[) 


xdetiVi" 1 , 


(9) 


where Zbath is the partition function of the bath. The 
matrix elements of M at (i, j) are given by the hy¬ 
bridization function A a ^ a . {r[ — Tj). 

We perform importance sampling with respect to the 
partition function using the weight 


w(d ai (T i), 

dr n 

:Tri oc 


npn 

d ai (n)dl(Tl) 


■■ ,d an (T„);d ai (r0,--- ,d a '(r')) 

e -/3Wioc Td an (T n )d f a , {r' n ) ■ ■ ■ 


detM 1 


( 10 ) 


The local trace is evaluated using the matrix formalism: 
e -rU\ oc ^ ^ i n Eq (IQ) are represented in the eigen- 
basis of Ho- [38] Note that we do not have to transform 
Ho to the new single-particle basis in terms of d and d). 


Effect of the off-diagonal hybridization functions 

We briefly discuss the effects of the off-diagonal hy¬ 
bridizations in the LDA+DMFT results. Figure [5] 
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compares the spin and orbital moments computed 
with/without off-diagonal elements at U = 2.5 eV. The 
behavior near the phase transition is affected by dropping 
the off-diagonal elements. In the calculation with diago¬ 
nal hybridization functions, the transition is weakly first 
order, as evidenced by the small hysteretic region. On 
the other hand, we observed a continuous transition in 
the simulation with off-diagonal hybridization functions. 


FIG. 5. (color online). Comparison of the spin and orbital 
moments computed with/without off-diagonal elements of the 
hybridization function. A hysteresis behavior is seen in the 
results computed with the diagonal hybridization function. 









